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ABSTRACT 

There are ubiquitous flows observed in the solar atmosphere of sub-Alfvenic 
speeds, however after flaring and coronal mass ejection events flows can become 
Alfvenic. In this Letter, we derive an expression for the standing kink mode 
frequency due to siphon flow in coronal loops, valid for both low and high speed 
regimes. It is found that siphon flow introduces a linear spatially dependent 
phase shift along coronal loops and asymmetric eigenfunctions. We demonstrate 
how this theory can be used to determine the kink and flow speed of oscillating 
coronal loops with reference to an observational case study. It is shown that the 
presence of siphon flow can cause the underestimation of magnetic field strength 
in coronal loops using the traditional seismological methods. 

Subject headings: magnetohydrodynamics (MHD) - - magnetic fields - - Sun: 
corona 



Introduction 



The measurement of plasma flow speed in the solar atmosphere has traditionally been 
estimated with Doppler shift using spectr ometers onboard, e.g., SOHO and Hinode (e.g., 



Brekke et al. 1997 



Winebarger et al.ll2002l ). Some attempts have also been made by track- 



ing features thou ght to be associated with flow using imagers (j Winebarger et al.l 12001 



Chae et al.l 120081 ). It has been found that sub-Alfvenic flows around of 100 km s 1 are 



ubiquitous in the corona. However, less frequent but faster flows have been detected in 
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the v icinity of flaring events and coronal mass ejections (e. g.. Ilnnes et al.ll2001t lHarra et al. 



20051 ). even into the Alfvenic regime of 10 km s (e.g., Ilnnes et al.ll2003l ). Some of the 



reported flows are of siphon type, i.e., the flow is unidirectional from one end of the loop 
to the other. Spec troscopic identifications of siphon flows can be found, for example, in 



Teriaca et all (120041 ) with SUMER/SOHO and iTian et all (120081 ) with SUMER/SOHO and 
EIS/HINO DE. Very clea r iden tifications of siphon flows through im aging observations were 
reported bv lDovle et aD tood ) with TRACE and ITian et aD J2009[ ) with STEREO. 



As well as generating fast flows, flares and CMEs can also cause standing kink oscillations 
in coronal loops (see e.g., lAschwanden et al.lll999t iNakariakov et al.lll999l ). so it is natural 
to develop magnetohydrodynamic (MHD) theory that models the interaction of these waves 
with flow. Since the dynamics of coronal plasma are dominated by magnetic fields, in 
general, flow direction will be magnetic field aligned. Using Doppler shift, there may be 
large uncertainty in estimating flow speeds along coronal structures such as loops due to the 
line of sight effects. In the most extreme case, if the direction of flow is perpendicular to the 
line of sight, there will be no flow detected at all. Furthermore, if a loop with field aligned 
flow is oscillating with a kink mode this adds further difficulty in measuring the flow speed 
using Doppler shift. In this Letter we want to demonstrate how the problem of estimating 
the flow speed in oscillating coronal loops can be addressed by implementing the technique 
of magnetoseismology. 

Coronal seismology is an indirect way to obtain the magnitudes of fundamental plasma 
parameters expl oiting the obse rved oscillatory propertie s of the solar atmosphere. This idea 
was proposed by lUchidal (J1970j) and lRoberts et al.l (11984) and up to now, has successfully pro- 
vided estimates of the magnetic field strength (Nakariakov fc OfmanlboOl ; Van Doorsselaere et al 



2Q0 8j) and density sc a le height (|Andries et al. 



2005 



Verth et al.ll2008l ) in the corona. Also, 



Arregui et al. ( 2007 ); Goossens et al. ( 20081 ) established lower limits to the sub- resolution 



transverse inhomogeneity in coronal loops, and upper limits to internal Alfven speeds. 

In this Letter, for the first time, we demonstrate how the flow speed along oscillating 
coronal loops can be determined from observations of standing kink modes. We theoretically 
investigate the effect of a constant siphon flow on stationary transverse waves and link this 
model to an observational case study. 



2. Equilibrium model 



Consider an equilibrium model of a cylindrical axis-symmetric flux tube of radius R 
with constant axial magnetic field B , and a density contrast of pi/p c . The subindex "i" 
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and "e" refer to the internal and external part of the tube, respectively. The length of the 
tube is L and it is assumed that inside there is a unidirectional constant flow U and no flow 
outside. 

The stationary solution is constructed by superposing two propagating waves traveling 
in opposite directions. The effect of flow on propagating waves causes a frequen c y Do ppler 



shift, a phenomena studied in the past by many authors, e.g., iGoossens et al.l ( 119921 ) . In 
particular it was found that, in the zero-/? limit, a valid approximation for the coronal 
plasma, frequency is given by the following expression, 



u Pi T T _l_ / A PiPe U 2 

u = k ■ U ± fcCkWl - t ; rx — 

Pi + Pe V (Pi + Pc) 2 c 



2 - 
k 



where Ck the kink speed for a static equilibrium within the thin tube r egime (R <C L) 



The corresponding eigenfunctions (see for example lEdwin fc Robertslll983l ) are expressed in 
terms of Bessel functions in the radial direction and have a sinusoidal dependence in the 
longitudinal direction. 

To find a stationary solution, we impose that the forward (+ sign) and backward ( — 
sign) waves must have the same frequency. This condition can be only satisfied if the 
corresponding longitudinal wavenumbers are different, and according to Eq. (JTJ are given by, 

k ± = " (2) 

-*-U±c k Jl - r J r k ^ 

Note that in general k + > 0, and k^ < (we do not consider KH-unstable modes here, i.e., 
we assume the square root is always positive). 

We concentrate on the z— dependence of the stationary solution, neglecting the radial 
dependence because we are assuming the thin tube limit. The study of the full MHD 
eigenvalue problem would require a more involved analysis. The transversal displacement of 
the tube axis is thus 

£(£, z) = A sin {ut — k + z) + B sin (ut — k_z) , (3) 

with A and B constants determined by boundary conditions. Imposing line-tying conditions 
at the footpoints, i.e., 



£(M = 0) = 0, 
&t,z = L) = Q, 



(4) 
(5) 
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from the first condition we find that A = — B. Using this relation and the second condition, 
we obtain the following dispersion relation 

1 Pi u 2 

J- — - i - ~~!Z 



where 



, v Pi+Pc c£ j 

u = k c k ===, (6) 

PiPc U 2 
(Pi+Pc) 2 T 



7T 

fc = n _ n = 1,2,.... (7) 



The natural eigenfrequencies of the flux tube are given by Eq. ([6]) (see iTarovanl 120091. for 



the eq uivalent result for a pure Alfven wave, and the recently published work of 



Ruderman 



(120101 ). for a more general study). This expression gives the explicit dependence of frequency 
on U. If we compare with the eigenfrequencies of the static case, u = /coCk, it turns out that 
flow always leads to a frequency reduction, i.e., an increase in the period of oscillation. 

We now turn our attention to the time-dependence of the amplitude of oscillations along 
the loop. Using the fact that A = —B and Eq. ([6]), we have that Eq. ([3]) can be written, 
after some algebra, in the following elementary form 

£ (t, z) = C sin k z cos (ut + kjjz) , (8) 

where C is an arbitrary constant, and we have introduced the wavenumber due to flow, 

ku = k ====. (9) 

pi + pcC k 1 _^iPcU^ 

V (Pi+Po) C k 

When there is no flow, U — 0, we recover the standing wave pattern of the static case, 
being a solution separable in time and space. When there is flow, the solution is no longer 
separable but the interpretation is straight forward. The term with the cosine in Eq. (jSJ) 
represents a propagating wave traveling in the opposite direction of the flow with an effec- 
tive wavenumber, which is simply linearly proportional to flow speed in the slow flow regime 
(U Ck). While the term with the sine is just the envelope of the traveling wave, satis- 
fying the line-tying conditions at the footpoints. Key to the present investigation, Eqs. jS]) 
and (Q show there are two main observational signatures of siphon flow in a coronal loop 
kink standing wave: 



Signature 1.- There is a linear phase dependence of the standing kink mode along the 

loop. 
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Signature 2.- In one full period the eigenfunctions of the standing kink mode will mostly 
exhibit an asymmetry about the center of the loop. 

If either Signature 1 or 2 are present in data, then this may be an indication there is a siphon 
flow present but if both are present this is a stronger indication of siphon flow. Figure [1] 
illustrates the dependence of amplitude along a coronal loop at different times due to constant 
siphon flow, illustrating Signature 2 of siphon flow. It is worth noting that in the absence of 
siphon flow Signature 1 may simply indicate the presence of a propagating kink wave. The 
correct interpretation of data must be done on a case by case basis, taking into account loop 
geometry and all estimated wave parameters. 



3. Expressions for siphon flow and kink speed 

Usually the data analysis of standing kink oscillations is based on a fit of the loop 
displacement at a given point along the loop of the following form 

f = £ cos(cut + (p)e~ t/TD . (10) 

The parameters that are determined from the fit are £q, u, <f), and td. The damping time, td, 
is not of interest in this work since a damping mechanism is most likely required to produce 
the attenuation of the signal with time. Regarding the other parameters it is important to 
point out that it is possible, depending on the data set, t o derive their values along some 



portion of the loop and not just at a single point (see e.g., IVerwichte et al.ll2010i ). 



Usually the fundamental mode is reported in the observations, mea ning that n = 1 



in Eq . (0) but there is also possible evidence of higher overtones (see e.g., IVerwichte et al. 



2004 ). The loop length, L, can be estimated using a circular shape or a 3D reconstruction, 



meaning that ko, is well determined (up to the uncertainties in L). 

When there is no flow, the determination of u (calculated using the fit) allows us, 
together with the value of ko, to calculate the kink speed of the tube, i.e., Ck = oj/ko. 
However, with only estimates for u and ko, we cannot determine the kink speed if flows are 
present. We also need additional information about the phase and plasma densities, p\ and 
p e . From the Signature 1 of siphon flow in Section [2J phase should show a linear dependence 
with position along the loop with gradient ky and this can be estimated with a linear fit to 
data. Also, the values of p x and p e can be estimated from observed intensity, since emission 
measure is approximately proportional to density squared if one considers the intensity of 
a clean emission line formed in the transition region or corona. In imaging observations, 
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usually a passband is wide enough to include several emission lines with different formation 
temperatures and the integrated intensity of this passband is no longer simply proportional 
to the density squared. The density diagnostics is possible if you have observations of density 
sensitive line pairs. But this usually needs spectroscopic observations (e.g., SUMER, CDS, 
EIS). 

Once u, k , kjj, p x and p e are determined from observation, it is straight forward to 
calculate the two magnitudes of interest, namely U and c\. We only need to combine Eq. ([6]) 
and Eq. © to find 

U = — -rrw, 11 



Pi ^0 



and 



k 2k% + Pl + Pc U + 2V kl \kl + ( Pi + Pe f , 



Note that Eq. (TTTT) contains the information about the flow direction as well, if ku > 
< 0) then we have that U > (U < 0) assuming that fco > 



4. Observational case study 



There have been a number of observational studies which have s hown a gradient in phas e 
along coronal loops osc i llatin g wit h the standing kink mode, e.g., IVerwichte et al.l (120041 ). 
Van Doorsselaere et al.l (120071 ) and lDe Moortel fc Brady! (120071 ). However, it is unclear from 
these studies whether the observed phase difference is due to siphon flow or wave propagation 
since there is only evidence of Signature 1 in these cases. For this reason we focus instead on 
the particular case study of IVerwichte et al.l (120101 ). since both Signature 1 and 2 are present 
in the data, providing a stronger argument for the si phon flow interpretation. Using combined 
TRACE and EIT observations, IVerwichte et al.l (120101 ) analyzed the standing kink mode 
generated in a coronal loop that was part of a large arcade. The coronal loop, positioned 
off-limb, starts oscillating after an Ml. 5 GOES level flare and associated CME occurs at 
about 11:04 UT. Regarding Signature 1, IVerwichte et al.l (120101 ) fitted a linear function to 
the observed phase along a 120 Mm portion of the loop (about 18% of the total loop length 
total L = 680 Mm) and this is shown in Figure 6 of their paper. Evidence of Signature 2, 
an asymmetry in the amplitude about the loop half length is shown in their Figure 8. 



Verwichte et al.l (120101 ) perform a fit to the loop displacement as in Eq. (fTOl) . The values 



of the parameters of interest (see Table I of their paper) are u = 2.60 x 10 3 s 1 (P = 2418 
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s), ko — 4.62 x 10 3 Mm (L = 680 Mm, and n — 1). The authors also provide a linear fit to 
the phase as a function of position along the loop (denoted by s) which is 



0.93 - 3.67 x 10 



-3 



IMm 



-80) [rad], 



(13) 



meaning that ki 



-3.67 x 10~ 3 Mm" 1 . Using Eq. (HI]) we find that U = -1445 km 



s -1 , while from Eq. ([12]) we obtain that Ck = 1610 km s _1 (assuming a density contrast of 
Pi/ p e = 5). Using the static model we find that Ck = 563 km s -1 , almost three times smaller 
than the inferred value assuming a siphon flow. The value of the flow speed is negative, 
m eaning that the flow wo uld be towards the footpoint shown in the lower image of Figure 1 
of IVerwichte et all ( boiol ). 



It is relatively simple to estimate the errors in the calculations of the flow and kink 
speeds. These errors are due to the uncertainties in the quantities that appear in Eqs. (fTTT) 
and ( TT2"|) . Typically the errors are 5L = 0.05L, 5P = 0.05P, 5(p c / p{) = 0.4 and also assume 
that 8ku = 0.10/cf/. Applying the error propagation formula in Eqs. (ITT]) and (TT2"]) we find 
that the flow is 1445±755 km s _1 , and the kink speed is 1610 ±903 km s _1 . In this example, 
the errors in the speeds are not small, and are specially sensitive to the uncertainty in k\j. It 
is worth noting that the calculated speeds and errors are based on the assumption of constant 
flow speed along the loop. It is clear that, for this particular example, the assumption of 
a constant siphon flow significantly increases the inferred kink speed (almost by a factor 
of 3). Therefore the presence of siphon flow has an importa nt implication for the metho d 
of estimating magnetic field strength along coronal loops by iNakariakov fc Ofmanl ( 120011 ) , 
which assumes a static equilibrium. In the present example the magnetic field strength would 
therefore be underestimated by a factor of 3. 

Assuming there is a siphon flow present, then the flow is in the fast (Alfvenic) flow regime 
of 10 3 km s _1 . The oscillation event analyzed here takes place in a coronal loop arcade in 
the vicinity of an M class flare and CME. In similar events, fast flow sig natures have often 



been measured in Dopp ler shift (e.g., Ilnnes et al.ll200ll ; lHarra et al.ll2005l ). After an X class 
flare Ilnnes et al.l (120031 ) measured similar Alfvenic flow speeds in the range 800 — 1000 km 



s _1 in a coronal arcade using combined SUMER and TRACE observations, and it must be 
emphasized these Doppler shifts are likely to be underestimates of true field aligned flows 
speeds due to line of sight effects. 

Alfvenic flow would cause significant asymmetry about the loop half length in the eigen- 
function (Signature 2) along the loop. Interestingly, the estimated eigenfunction derived by 
Verwichte et al.l (120101 ) using EIT, shows that the fundamental mode has s uch an asymme- 



try. T his is clearly seen in Figure 8(a) (see the region s/L G [0,0.4]) where IVerwichte et al. 



( 120101 ) fit the observed displacement (solid line) with that predicted assuming a planar loop 
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model (dashed line). However, according to our model, this asymmetry can be explained by 
fast axial flow (see Figured]). 



5. Summary and Conclusions 



Using a simple model of a coronal loop, we have derived an analytical expression for the 
frequency of oscillation of a thin tube with a constant siphon flow, valid for both slow and 
fast flow regimes. It was found that calculation of the kink speed and estimation of magnetic 
field assuming a static model give values which are smaller than the ones obtained in the 
presence of flow. This questions the validity of traditional seismological methods assuming 
a static background to estimate the magnetic field along coronal loops. 

Contrary to the static case, different posi tions along the tu be oscillate with a different 
phase, in agreement with the recent results of iRudermanl (120101 ) . However, we have pointed 
out in our work the importance of the phase shift, since it has a linear behavior with the axial 
coordinate along coronal loops and in the low flow regime is simply linearly proportional to 
the flow speed. Furthermore, the presence of flow breaks the symmetry of the eigenmodes 
about the loop half length. 

According to our m odel, assuming the line ar phase shift reported along the post-flare 
coronal loop analyzed by IVerwichte et al.l (120101 ) is caused by a siphon flow, the flow would 
be in the Alfvenic regime. This is not completely unreasonable since, in the dynamic post- 
flare environment of coronal loops s uch fast flows have been estimated using Doppler shift 
and feature tracking techniques (e.g. Jlnnes et al.ll200ll . 120031 ; lHarra et al.ll2005l ). In any case, 
other cases need to be investigated in order to prove that this technique can be also used 
with sub- Alfvenic flows. 

The most important observed parameter that determines the value of the flow (and kink 
speed) in our model is the value of ky. The value of this parameter was estimated using 
only a small portion of the loop. If we determine the phase dependence along the whole loop 
(SDO should help in this regard) we would have a much more accurate determination of the 
speeds and a significant reduction in the errors. This more complete information will also 
allow us to test whether the flow is unidirectional or not, according to the sign of kjj at the 
footpoints. 

In reality flows along coronal loops are much more complex than the simple model 
proposed in this Letter. In particular, they are probably non-symmetric r espect to the loop 



apex, and more important, they are usually transient (see for example, iDoyle et al.l 12006 



Tian et al.l 120091 ) . Nevertheless, we must understand the simplest possible configuration 
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before additional complications are added. In this regard, we have made the first step in 
assessing the effect of constant siphon flows on standing kink oscillations. 
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Fig. 1. — Time evolution of the fundamental mode (n = 1) as a function of position along the 
loop for the case with U = 0.5ck- The curves represent the eigenfunction at different times (the 
interval is P/8, where P = 2tt/u). The continuous lines represent motion of the loop in the negative 
direction, while dashed lines represent motion in the positive direction. 



